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Two dimensional materials are important for electronics applications. A natural way for electronic 
structure engineering on two dimensional systems is on-plane chemical functionalization. Based on 
density functional theory, we study the electronic structures of fluorine substituted planar polysilane 
and graphane. We find that carbon and silicon present very different surface chemistry. The indirect 
energy gap of planar polysilane turns to be direct upon fluorine decoration, and the gap width is 
mainly determined by fluorine coverage regardless of its distribution on the surface. However, 
electronic structure of fluorine doped graphane is very sensitive to the doping configuration, due to 
the competition between antibonding states and nearly-free-electron (NFE) states. With specific 
fluorine distribution pattern, zero-dimensional and one-dimensional NFE states can be obtained. 
We have also studied the chemical modification with -OH or -NH2 group. Carbon seems to be too 
small to accommodate big functional groups on its graphane skeleton with a high concentration. 



I. INTRODUCTION 

Low-dimensional materials are expected to play a key 
role in nanotechnology. For example, it is desirable to use 
two dimensional materials to build electronics on flexible, 
lightweight, and cheap substrates like plastics.^ Silicon is 
the fundamental semiconducting material in electronics. 
Due to its relatively large size, it prefers to form four co- 
valent a bonds in tetragonal coordination (sp^ hybridiza- 
tion) instead of forming strong planar tt bonds through 
sp^ hybridization. A stable two dimensional silicon mate- 
rial is planar polysilane,^ which is composed of a buckled 
hexagonal silicon framework and hydrogen atoms on both 
sides to saturate dangling bonds. Layered polysilane has 
been synthesized experimentally^^^ and due to confine- 
ment effect, its electronic structure is different from bulk 
silicon 

After the successful experimental preparation of 
graphene^ carbon based materials become a rising star 
for two dimensional electronics.^^ Compared to sili- 
con, carbon have smaller radius, and stable tt bonds 
can be formed between carbon atoms within a plane. 
Therefore, the stablest two dimensional carbon struc- 
ture is a hexagonal lattice, namely graphene. Carri- 
ers in graphene show massless Dirac Fermion behavior, 
and there is a symmetry between electron and hole in 
graphene. Since graphene is a semimetal, an energy gap 
must be opened before any 'graphenium' microproces- 
sor can be realized. A natural way for this purpose is 
to convert the sp^ bond in graphene to sp^ bond via 
hydrogenation. Such a carbon based sp^ structure is 
called graphane, the carbon analogue of planar polysi- 
lane, which was theoretically predicted to have a finite 
energy gap.^^^^^ Recently, graphane has been synthesized 
by hydrogen plasma treatment or by electron-induced 
dissociation of hydrogen silsesquioxane.^^ 

For various electronics applications, it is very de- 
sirable to have the capability to fully control the en- 
ergy gaps and electronic structures of planar polysi- 
lane and graphane. On-plane chemical functionaliza- 



tion is an important way towards this direction. Hy- 
droxy group^iiii^ substitution on planar polysilane has 
been studied. The product, siloxene, has a direct energy 
gap instead of the indirect one in unsubstituted planar 
polysilane i^ii^ For two dimensional carbon based mate- 
rials, previous electronic structure engineering attempts 
by chemical means mainly focuses on graphene nanorib- 
bon edge modification^i^iiiii^ However, on-plane func- 
tionalization does not need to cut the graphene sheet, 
and it is expected to be a more feasible way to realize 
electronics on graphene (for example, by printing circuit 
on graphene). Although graphene oxide (GO) represents 
an on-plane chemically modified carbon material, it is 
mainly studied as an intermediate for massively produc- 
tion of graphene and the detailed structure of GO 
is still unclear.?^^^^^ 

Despite its importance, a systematical study of on- 
plane chemical functionalization for two dimensional sil- 
icon and carbon materials is still unavailable. Most im- 
portantly, how the electronic structure depends on the 
adsorption configuration of functional groups has not 
been discussed, which is very important for practical real- 
ization of electronic structure engineering in experiment. 
Another not considered important issue is the nearly- 
free-electron (NFE) states. NFE states are delocalized 
on the surface outside the atomic centers. When they 
are close to the band edge, NFE states become a deter- 
mining factor of the band gap width j^^i^^i^^ Since NFE 
states have unique space distribution, we expect their 
response to surface chemical functionalization to be dif- 
ferent than other states. 

In this article, we systematically studied the on-plane 
chemical modification on planar polysilane and graphane. 
Functionalization leads to significant electronic structure 
modification. Interestingly, electronic structure of fluori- 
nated polysilane is much less sensitive to the F adsorp- 
tion configuration than graphane. The remainder of the 
paper is organized as the following. In section II, we de- 
scribe the computational details. The main results are 
discussed in section HI. We first compare the electronic 
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structure of two dimensional polysilane and graphane, 
and the NFE state on graphane is introduced. Then, 
we discuss different behaviors upon F substitution for 
these two systems. Other functional groups are finally 
discussed. Section IV concludes the paper. 



II. COMPUTATIONAL DETAILS 

To investigate the geometric structures and elec- 
tronic states of polysilane and graphane based systems, 
we performed first-principles DFT calculations with 
the Vienna Ab initio Simulation Package (VASP);^i^ 
Projector-augmented wave (PAW) method^^^^^ was used 
for electron- ion interaction, and the Perdew-Wang forn>2^ 
of the generalized gradient approximation (GGA) was 
adopted for electronic exchange and correlation. A 
kinetic-energy cutoff of 400 eV was selected for plane 
wave basis set. For pristine planar polysilane and 
graphane we chose a 1x1 unit cell, and a 2x2 supercell 
was used to consider chemical modification. The vacuum 
space between two neighboring two-dimensional sheets is 
larger than 15 A, and dipole correction is applied to com- 
pensate artificial dipole interaction between neighoring 
supercells in the direction perpendicular to the plane. A 
13x13x1 /c-point mesh was used in geometry optimiza- 
tions, and a 27x27x1 grid was used for static electronic 
structure calculations. Geometry structures were relaxed 
until the force on each atom was less than 0.01 eV/A 
and the the convergence criteria for energy is 10~^ eV. 
A series of unit cells were scanned for each system to 
obtained the lattice parameters with the lowest energy. 
Charge populations were calculated using Bader's atom 
in molecule (AIM) method based on charge density topo- 
logical analysis 




FIG. 1: Schematic (a) side view and (b) top view of planar 
polysilane and graphane. Silicon or carbon is in red, and 
hydrogen in blue. An 2x2 unit cell is marked in the top view. 



two dimensional polysilane, with the valance band max- 
imum (VBM) at the F point and the conduction band 
minimum (CBM) at the M point. The direct band gap 
at the F point is just 0.12 eV broader. Graphane has an 
about 3.43 eV direct gap at F, as also suggested by Sofo 
et al^ 
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FIG. 2: (a) Band structure of planar polysilane, and corre- 
sponding charge densities of (b) the VBM and (c) the con- 
duction band edge at the F point, (d) Band structure of 
graphane, and corresponding charge densities of (e) the VBM 
and (f) the CBM at the F point. In the band structures, oc- 
cupied bands are in black, and unoccupied bands are in blue. 



III. RESULTS AND DISCUSSION 

A. Geometrical and electronic structures of planar 
polysilane and graphane 

First, we compare the geometrical and electronic struc- 
tures of planar polysilane and graphane. Both of them 
are composed of a buckled hexagonal plane of silicon or 
carbon, and hydrogen atoms on both sides of the plane 
to saturate the sp^ a bonds (Fig. [1]). The optimized Si- 
Si bond length is 2.36 A, similar to its calculated bulk 
value (2.35 A). Si-H bond length is 1.50 A. The opti- 
mized C-C and C-H bond length in graphane is 1.54 and 
1.11 A, respectively. The C-C bond length in graphane 
is almost identical to its value in bulk diamond, and it 
is also in agreement with the previous theoretical results 
on graphane. ^^'^^ 

Because its silicon or carbon atoms are sp^ hybridized, 
both planar polysilane and grphane have a significant en- 
ergy gap (Figure [2]). In agreement with Hirayama et al.,^ 
we obtained a quasi-direct band gap about 2.20 eV for 



The VBM of polysilane is a Si-Si bonding state. Due 
to symmetry, it is two- fold degenerated at the F point. 
The conduction band edge at F point is the antibond- 
ing state between Si 3pz and H Is orbitals, mainly dis- 
tributed on Si. VBM of graphane have similar bond- 
ing character with planar polysilane. However, CBM 
state of graphane is notably different from that of ploysi- 
lane, where a strongly delocalized feature exists above 
and blow the graphane plane. Such a delocalized state 
is known as NFE state*^ When it is pushed to the band 
edge, NFE state becomes a determining factor of the 
width of band gap. We note that the NFE-like char- 
acter of its CBM state has not been noticed in previous 
study on graphane. 

Why CBM of graphane is an NFE state, while that 
of polysilane is not? One reason is that the antibond- 
ing states of graphane are higher in energy than those 
of planar polysilane. On the other hand, NFE states of 
graphane are also strongly stabilized compared to those 
of planar polysilane. In fact, at the F point, the lowest 
NFE states of polysilane and graphane locate at about 
0.7 and 1.4 eV below the vaccum level, respectively. 
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There are two possible reasons for the stabilization of the 
NFE states in graphane. First, compared to notably neg- 
atively charged hydrogen in polysilane, hydrogen atoms 
in graphane are almost neutral. A negatively charged 
surface is repulsive to electrons, and thus upshift NFE 
states. Second, the radius of carbon is smaller than sil- 
icon, which leads to a much smoother graphane surface 
than the polysilane surface. Actually, the distance be- 
tween two hydrogen atoms in graphane is 1.35 A shorter 
than that in polysilane. A smoother surface is expected 
to stabilize NFE states. 



B. Fluorinated planar polysilane 

Fluorine has a much larger electronegativity than hy- 
drogen. Therefore, the substitution of hydrogen in pla- 
nar polysilane by fluorine may lead to a significant elec- 
tronic structure modification. We first consider the sim- 
plest case where all hydrogen atoms are replaced by fiu- 
orine atoms (100% F coverage). F substitution does not 
change the Si-Si bond length, and the optimized Si-F 
bond length is 1.63 A. 

The band structure of F-substituted polysilane (Fig. 
[3li) is very different from the unsubstituted system. The 
band gap (0.62 eV) becomes much smaller, and it is a di- 
rect gap at the F point. Therefore, F-substituted polysi- 
lane will have a much higher luminous efficiency than 
bulk silicon. CBM of F-substituted polysilane at the F 
point is mainly an antibonding state between Si and F 
molecular orbitals. Similar to siloxene,^ the indirect gap 
to direct gap transition from unsubstituted to substituted 
polysilane may caused by stronger mixing with F 25 state 
at r point, which pushes down the CBM at F to lower 
than the local minimum at the M point. VBM is still 
mainly a Si- Si bonding state. However, compared to pla- 
nar polysilane (Figure [2]3), there is a small hybridization 
with the F 2s orbital. 




FIG. 3: Band structure of fluorinated planar polysilane for the 
most stable configuration with a (a) 25%, (b) 50%, (c)75%, 
and (d) 100% F coverage, respectively. 



It is also interesting to check the electronic structure 
evolution with the F coverage. Using a 2x2 supercell, 
we consider 25%, 50%, and, 75% F substitutions, which 
corresponds to 2, 4, and 6 ffuorine atoms in the supercell, 
respectively. For 25% and 75% F coverages, we consider 
all possible F substitution configurations. For 50% F 
coverage, 6 configurations are considered. All calculated 
configurations for different F concentrations are plotted 
in Fig. m 




FIG. 4: Fluorine substitution configurations. Fluorine atoms 
are in red, and hydrogen in blue. Black lines represent the 
silicon or carbon skeleton. 

The energetically most stable configuration with dif- 
ferent F coverage is 25-c, 50-c, and 75-c, respectively. 
However, the energy difference between different config- 
urations for a certain F concentration is very small, typ- 
ically within 0.1 eV. At the same time, with a certain F 
coverage, configuration difference has little effect on the 
band structure. As shown in Fig. [31 band gap decreases 
with the increase of the F coverage (Table [Ij). 

In all configurations, CBM is mainly of antibonding 
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character between Si-F or Si-H molecular orbitals. Since 
the Si-F antibonding state is lower in energy than the 
Si-H antibonding state (refer to the 100% F coverage 
case and the unsubstituted case), more weight of Si-F 
antibonding character leads to lower energy of the band. 
Therefore, as the percent of F substitutions increases, the 
energy of CBM significantly decreases. 

VBM of fiuorinated polysilane is a Si-Si bonding state 
with a small hybridization with F orbitals. The energy 
of VBM also slightly decrease with the F coverage, which 
can be understand by charge transfer between Si and F. 
Since fluorine has a larger electronegativity than silicon, 
an electron-transfer from Si to F is expected. Charge 
population Bader analysis shows that the amount of elec- 
trons transfer from Si to F increases with the increase of 
the F concentration (Table [11). As VBM is mainly a Si-Si 
bonding state, a positively charged Si skeleton will push 
down the VBM. Because the downshift of CBM is more 
significant than VBM, we still get smaller bang gap for 
strongly fiuorinated polysilane. 



C. Fluorine functionalization on graphane 

If all hydrogen atoms in graphane are substituted by 
fiuorine, the C-C bond length will be slightly increased 
by 0.03 A to 1.57 A, and the optimized C-F bond length 
is 1.38 A. As shown in Fig. O F substitution decrease the 
band gap of graphane to about 3.10 eV. CBM at the F 
point is an antibonding state between C and F orbitals, 
and VBM is a mix of C-C bonding state and F 2s or- 
bital. Compared to VBM of fiuorinated two dimensional 
polysilane, the hybridization with the F orbitals here is 
stronger. 
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FIG. 5: (a) Band structure and charge densities of the (b) 
CBM and (c) VBM at F for fiuorinated graphane. 

We also study different F coverage with all configura- 
tions listed in Fig. [H For silicon, when the F coverage is 
fixed, different substitution configurations lead to simi- 
lar band structures. Interestingly, the situation is totally 
different for F partially substituted graphane, where the 



F ordering plays a key role in the electronic structure of 
the doped systems. 

With the 25% F coverage, we have three configura- 
tions. The 25-c configuration has similar band structure 
with 25-b. Therefore, in Fig. [6l we only plot the band 
structures for 25-a and 25-b. There is an about 1.5 eV dif- 
ference between the energy gaps of 25-a and 25-b. VBM 
in both configurations is mainly a C-C bonding state with 
a small mixing of F 2p orbital. Just like F-substituted 
planar polysilane, its position slightly decreases with the 
F coverage. The CBM of 25-a is significantly lower than 
that of 25-b. CBM is mainly a NFE state similar to that 
of unsubstituted graphane (Fig. [2^). Since the NFE state 
is delocalized apart from the system, its energy can be 
easily affected by a shift of local vacuum level. A dipole 
layer on surface can modify the local vacuum level. In 
Table m we list the dipole moment for different configu- 
rations. Since, in 25-b and 25-c, the two fiuorine atoms 
is distributed on two sides of the graphane plane, there 
is no net dipole moment for these two configurations. 
However, in 25-a, the two fiuorine atoms is on the same 
side. Because fiuorine atoms get electrons from the car- 
bon skeleton, there is a dipole moment perpendicular to 
the surface, which lower the local vacuum level and thus 
also the NFE conduction band. 




K TMKK TMKK TMKK TMKK TMK 
50-a 50-b 50-c 50-d 50-e 



FIG. 6: Band structures of fiuorinated graphane with different 
configurations. 

For 50% coverage, 50-f has similar electronic structure 
with 50-e, and will not be discussed here. In the 50-a 
configuration, all F atoms are on the same side, and the 
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TABLE I: Energy gap(£^p), averaged charge population (p) on Si or C, and dipole moment (/i) of planar polysilane and graphane 
with different degree of fluorine functionalization. Eg is in eV, p in electron, and /i in ex A per supercell. 










25-a 


25-b 25-c 


50-a 


50-b 50-c 


50-d 50-e 50-f 75-a 


75-b 75-c 100 


Si-F 


Eg 


2.20 


1.94 


1.86 


1.96 


1.67 


1.37 


1.62 


1.52 


1.49 1.57 


1.13 


1.06 


1.21 0.62 




P 


3.46 


3.37 


3.38 


3.38 


3.32 


3.32 


3.32 


3.32 


3.32 3.32 


3.26 


3.27 


3.26 3.19 









0.30 








0.60 











0.30 0.30 


0.30 








C-F 


Eg 


3.43 


3.43 


4.94 


4.86 


2.60 


4.56 


5.44 


5.16 


4.54 4.46 


4.17 


4.14 


4.77 3.10 




P 


3.98 


3.87 


3.88 


3.88 


3.71 


3.73 


3.74 


3.74 


3.73 3.73 


3.57 


3.58 


3.59 3.41 




P 





0.36 








0.64 











0.32 0.32 


0.28 









NFE-like CBM (distributed on the H side) is lowered by 
the dipole layer generated by charge transfer between C 
and F. As a result, the 50-a configuration has the smallest 
energy gap in all 50% F coverage configurations. 

CBM of the 50-b configuration is an antibonding state 
between C and F molecular orbitals and the next two 
bands higher in energy at the F point are mainly NFE 
states (Fig. [7^ and[7)3). Here, the lowest NFE state is 
higher than the lowest antibonding state. This is because 
that the F doping destroy the smooth H surface, which 
destabilize the NFE state. There is no dipole moment 
in 50-b, and the dipole stabilization mechanism does not 
help here. Interestingly, due to the chain configuration 
of H atoms, the NFE state forms one-dimensional distri- 
bution. One dimensional NFE state has potential appli- 
cation in electronics as an ideal transport channel when 
it is close to the Fermi level. 



(a) 


(c) 


m 


(b) 


(d) 






ft ft 





FIG. 7: (a) Top view and (b) side view of the lowest NFE 
state of 50-b. (c) Top view and (d) side view of the lowest 
NFE CBM state of 50-e. 



Surprisingly, an indirect band gap is obtained for 50- 
c. At the r point, the lowest unoccupied band (LUB) is 
mainly a ribbon-like NFE state, but also with a strong 
hybridization with some C-H and C-F anti-bonding char- 
acter. While the CBM at the /c-point between F and M 
is mainly of antibonding character. This fact indicates 
a strong interaction between the NFE state and the an- 
tibonding state, which may be a reason for the indirect 
gap. In 50-d configuration, we also have hydrogen chains 
formed, which leads to two one-dimensional NFE bands 



(LUB+2 and LUB+3). However, in this configuration, 
hydrogen chains on the two sides of the graphane plane 
are not parallel. 

Due to the dipole moment effect, CBM of 50-e is also a 
NFE state. Because neighboring H atoms form triangles, 
NFE states in 50-e are also of the shape of triangle (Fig. 
[Tt and[71i). Such zero-dimensional NFE states demon- 
strate the flexibility to tune the distribution and energy 
of NFE states over two dimensional surface by chemi- 
cal functionalization. Due to the smaller dipole moment 
compared to 50-a, the NFE state in 50-e is higher, which 
leads to a larger gap. 

The band structures of 75-a and 75-b are shown in Fig. 
[6l and the band structure of 75-c is similar to 75-b. At 
the 75% F coverage, due to the low coverage of hydro- 
gen, NFE state is not very stable, and CBM for all the 
three configurations is antibonding states between C and 
F molecular orbitals. Since the antibonding state is less 
affected by dipole moment, the gap difference between 
different configurations at 75% F coverage is much less 
than those at 25% or 50% F coverages. 



D. Other functional groups 

Hydroxy (-0H) and amine (-NH2) groups, both of 
which have multiple atoms, are also considered. For sim- 
plicity, we only discuss their 100% coverage on both sides 
and 50% coverage on one side (full coverage on that side, 
the 50-a configuration). 

OH or NH2 functionalization does not change the lat- 
tice parameters of two dimensional polysilane, which is 
consistent with the relatively large Si- Si bond length 
and the week steric hindrance between functional groups. 
However, for smaller carbon, due to the repulsion be- 
tween substituted radicals, the optimized C-C bond 
length upon OH or NH2 substitution is notably larger 
than that of graphane. With 50% and 100% OH cover- 
age, the C-C bond length is 1.57 and 1.62 A, respectively. 
For NH2 functionalization, the corresponding C-C bond 
length is 1.62 and 1.68 A, respectively. The elongated C- 
C bond length indicates that it is difficult for graphane to 
accommodate large functional groups at a relatively high 
coverage, as also suggested by a previous GO structure 
study. 
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TABLE II: Energy gap(^^), averaged charge population (p) 
on Si or C, and dipole moment (/i) of planar polysilane and 
graphane with different degree of OH and NH2 functionaliza- 
tion. Eg is in eV, p in electron, and /i in ex A per supercell. 
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As shown in Fig. [8l the energy gap is also decrease with 
the increase of the functional-group coverage for both OH 
and NH2 functionalized polysilane (Table However, 
contrast with the F substitution case, VBM rise with the 
increase of the functional-group coverage. This may due 
to stronger mixing between Si-Si bonding states and N or 
O 2p orbit als, since they are higher than F 2p orbital and 
more close to the Si 3p orbital. The strong hybridization 
between Si and O or N orbitals is consistent with the 
split VBM, which is degenerated in fluorinated polysi- 
lane. CBM at the F point is mainly of the antibonding 
character between Si-0 or Si-N molecular orbitals. 




FIG. 8: Band structure of (a) 50% and (b) 100% OH covered 
or (c) 50% and (d) 100% NH2 covered planar polysilane. 

For both OH and NH2 doped graphane, we obtain a 
smaller gap for the 50% coverage compared to the 100% 
case. This is because the NFE state on the hydrogen 
atom in the 50% coverage case is lowered by the dipole 
layer developed in the system. In all four cases, VBM 
is a C-C bonding state mixed with O or N 2p orbitals. 
The mixing is very strong, which leads to an even larger 
splitting of VBM at the F point compared to OH or NH2 
substituted planar polysilane. The position of VBM is al- 
most not affected by the functional group coverage, which 
indicate that the N and O 2p orbitals are close in energy 
to the C-C bonding state. 




FIG. 9: Band structure of (a) 50% and (b) 100% OH covered 
or (c) 50% and (d) 100% NH2 covered graphane. 



CBM of 100% OH covered graphane is also the an- 
tibonding states of C and O molecular orbitals. But 
for 100% NH2 covered graphane, CBM has a significant 
NFE-like character due to the hydrogen atoms in the rad- 
icals (Fig. [TQ|) . Such kind of NFE states provide more 
flexibility to use NFE state to tune the electronic and 
optical properties of two-dimensional materials. 




FIG. 10: CBM density in 100% NH2 covered graphane. 



IV. CONCLUSIONS 

In summary, we have calculated the on-plane chemical 
modification of two dimensional polysilane and graphane. 
Chemical functionalization can change the indirect gap 
in planar polysilane to a direct gap, which will strongly 
modify the optical properties of the system. For fluori- 
nated polysilane, the gap decreases with the F coverage, 
independent on the F adsorption configurations. How- 
ever, the electronic structure of fluorinated graphane is 
sensitively dependent on the F doping geometry. This 
is because the NFE states of graphane is close to the 
conduction band edge. Energy of NFE states is more 
sensitive to the local vacuum potential and the doping 
structure. Functionalization by OH and NH2 groups is 
also discussed. 
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